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Abstract. We prove stability estimates for the ENO reconstruction and ENO interpolation proce- 
dures. In particular, we show that the jump of the reconstructed ENO pointvalues at each cell interface 
has the same sign as the jump of the underlying cell averages across that interface. We also prove 
that the jump of the reconstructed values can be upper-bounded in terms of the jump of the underly- 
ing cell averages. Similar sign properties hold for the ENO interpolation procedure. These estimates, 
which are shown to hold for ENO reconstruction and interpolation of arbitrary order of accuracy and on 
non-uniform meshes, indicate a remarkable rigidity of the piecewise-polynomial ENO procedure. 
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1. Introduction and statement of main results 

The acronym ENO in the title of this paper stands for "Essentially Non-Oscillatory", and it refers 
to a reconstruction procedure, which generates a piecewise polynomial approximation of a function 
from a given set of its cell averages. The essence of the ENO procedure, which was introduced by 
Harten et. al. in [[1211 . is its ability to accurately recover discontinuous functions. The starting point is 
a collection of cell averages {v,},^^ over consecutive intervals // = 1/2,-^/+ 1/2)' 

(1.1) Vi := 777 / v{x)dx, 

I A- 1 Jii 
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from which one can form the piecewise constant approximation of the underlying function v{x), 

But the averaging operator £/v{x) is limited to first order accuracy, whether v is smooth or not; for 
example, if v has bounded variation then ||v — £/v||^i = ^(/z). The purpose of the ENO procedure 
(abbreviated by M) is to reconstruct a higher order approximation of v{x) from its given cell averages, 

(1.2) ENO: s^vix) = J^v^l/, W ^ Wv(x) := Y^fkWM- 

k k 

Here, fk{x) are polynomials of degree p — I such that the piecewise-polynomial ENO reconstruction 
^£/v{x) satisfies the following two essential properties. 

Accuracy: First, it is an approximation of v{x) of order p in the sense that 

(1.3) Mj^v{x) = v{x) + ff{hP), 

where h = max,- Typically, the requirement for accuracy is sought whenever v(-) is suffi- 
ciently smooth in a neighborhood of x. Here, however, ( |1.3 ) is also sought at isolated points 



of jump discontinuities. Thus, if we let v(jc,^i/2+) and v(x,_|_i/2— ) denote the point-values of 



v{x) at the left and right of the interface at jc,_|_i/2, then ( 1.3 ) requires that the corresponding re- 
constructed point- values, vT^^^ := M£/v{xi^i/2—) = fi{xi+i/2) and v^^^ •= ^■^v{xi^i/2+) = 
fi+ii^i+i/i)^ satisfy 

l^r+i/2-^(-^'+V2-)l + l^ii/2-^(^^-+V2+)l = ^(^'')- 
To address this requirement of accuracy, the /,'s are constructed from neighboring cell av- 
erages {v,:+j}y^^^^ for some k G {—p -|- 1, ... ,0}. The key point is to choose an adaptive 
stencil, 

i ^ {vi+k: ■ ■ ■ , V/, ■ ■ ■ , Vi+k+p^i}, 
based on a data-dependent shift k = k{i). This enables the essential non-oscillatory property 



( 1.3 1, while making the ENO procedure essentially nonlinear. 



Conservation: The second property sought in the ENO reconstruction is that the piecewise-polynomial 
ENO approximation be conservative, in sense of conserving the original cell averages, 

(1.4) -r^ f ^£/v{x)dx = Vi. 

The conservative property enables us to recast the ENO procedure in an equivalent formula- 
tion of nonlinear interpolation. To this end, let V{x) := / v{s)ds denote the primitive of 

J — oo 

v{x). The given cell averages {v/} now give rise to a set of point- values {V^+i/jjjeZ, 

(1.5) V.^y^-= I v{s)ds= l '^'^v{s)ds= |4|v^. 

A second-order approximation of these point-values is given by the piecewise linear inter- 
polant ^V{x) := Zk |^ (Vfc- 1/2(^^+1/2 +^fc+i/2(-«--«yt-i/2)) 1/* W- The ENO approxima- 
tion, Y.kP'k{x)hi,ix), is a higher-order accurate piecewise-polynomial interpolant, 

(1.6) ENO: ^V(x) ^ M^V{x) ■.= YPk{x)\{x). 
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It interpolates the given data at the nodes, i^(jc,±i/2) = ^±1/2' ^^d it recovers V{x) to high- 
order accuracy at the interior of the cells, ^^V{x) = V{x) + &{hP^^). Now, let Fi{x) be the 
unique p-th order polynomial interpolating the p+\ pointvalues . . . , V,+r+p for some 



shift r which is yet to be determined. Then, it is a simple consequence of ( 1.5 1 that / satisfies 



( [13] ) and ( [L4| ) and that 

F,i^) = V,-i/2 + fi{s)ds. 



-1/2 

In this manner, ENO reconstruction of cell averages is equivalent to ENO interpolation of the 
pointvalues of its primitive. We shall travel back and forth between these two ENO formula- 
tions. 

1.1. ENO reconstruction. When the underlying data is sufficiently smooth, the accuracy require- 
ments can be met by interpolating the primitive V on any set of ;? + 1 point- values 

Here, r is the (left) ojfset of the interpolation stencil, which is indexed at half-integers, to match the 



cell interfaces. To satisfy the conservation property (1.4), the stencil of interpolation must include 
V,_i/2 and V,^i/2. There are p such stencils, ranging from the leftmost stencil corresponding to an 
offset r = —p + 1/2 to the rightmost stencil corresponding to an offset of r = —1/2. Since we are 
interested in approximation of piecewise smooth functions, we need to choose a carefully shifted 
stencil, in order to avoid spurious oscillations. The main idea behind the ENO procedure is the use 
of a stencil with a data-dependent offset, r — r(i), which is adapted to the smoothness of the data. 
The choice of ENO stencil is accomplished in an iterative manner, based on divided differences of the 
data. 

Algorithm 1.1 (ENO reconstruction algorithm: selection of ENO stencil). 

Let point values of the primitive V;_p_|_i/2, • ■ • , Vi+p_i/2 be given, e.g., ( |1.5[ ). 

• Setri = -I/2. 

• For each j = 1, — 1, do: 

fif \V [xi+rj-U- ■ ■ ^Xi+rj+j] I < \V [Xi+rj,- ■ ■ .Xi+rj+j+l] \ ^ SCt rj+i =rj-\, 

1 Otherwise t-)- set r^+i = r^. 

• Set Fi{x) as the interpolant ofV over the stencil {VJ+yt}^^^. 

• Compute fi{x) :=F-{x). 

The divided differences V[xk,.. . are a good measure of the jth order of smoothness of V{x). 

Thus, the ENO procedure is based on data-dependent stencils which are chosen in the direction of 
smoothness, in the sense of preferring the smallest divided differences. 

The ENO reconstruction procedure was introduced in 1987 by Harten et. al. [[T2ll in the context of 
accurate simulations for piecewise smooth solutions of nonlinear conservation laws. Since then, the 
ENO procedure and its extensions, [|23l[T3ll24l[T4l[T5l[T6l . have been used with a considerable success 
in Computational Fluid Dynamics; we refer to the review article of Shu [25 1 and the references therein. 
Moreover, ENO and its various extensions, in particular, with subcell resolution scheme (ENO-SR), 
[[T3l , have been applied to problems in data compression and image processing in [[TTl [H [6l [211 111 III 
[21 [31 and references therein. 

There are only a few rigorous results about the global accuracy of the ENO procedure. In [[3, 
the authors proved the second-order accuracy of ENO-SR reconstruction of piecewise- smooth 
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data. Multi-dimensional global accuracy results for the so-called ENO-EA method were obtained in 
[|3i . Despite the extensive literature on the construction and implementation of ENO method and its 
variants for the last 25 years, we are not aware of any global, mesh independent, stability results. This 
brings us to the main result of this paper, stating the stability of the ENO reconstruction procedure in 
terms of the following sign property. 

Theorem 1.1 (The sign property). Fix an integer p> I. Given the cell averages {v,}, let^%£i^v{x) be 
the p-th order ENO reconstruction of these averages, as outlined in Algorithm \\ .l\ 

Wv(x) =£A(x)l4W, degAW <p-l. 



Let v^y^ '■— ■^■«2/v(X;_|_y2+) and v.^y^ '■— ^■^^{^i+i/2~) denote left and right reconstructed point- 
values at the cell interface JC/+i/2- Then the following sign property holds at all interfaces: 

if v,+i— v,>0 then 
if v,+i— v/<0 then 



(1.7) 



^■+1/2- 



^■+1/2 



In particular, if Vi^\ = then the ENO reconstruction is continuous across the interface, v. 



+ 1/2 

Moreover, there is a constant Cp, depending only on p and on the mesh-ratio |/|/,| of 



^■+1/2- 

neighboring grid cells, such that 



(1.8) 



0< 



1+1/2 ''i+l/2 
Vi+l -Vi 



The sign property tells us that at each cell interface, the jump of the reconstructed ENO pointval- 
ues cannot have an opposite sign to the jump in the underlying cell averages. The sign property is 
illustrated in Figure [T] which shows a third-, fourth- and fifth-order ENO reconstruction of randomly 
chosen cell averages. Even though the reconstructed polynomial may have large variations within 
each cell, its jumps at cell interfaces always have the same sign as the jumps of the cell averages. 
Moreover, the relative size of these jumps is uniformly bounded. We remark that the inequality on 
the left-hand side of ( |1.8| ) is a direct consequence of the sign property ( |1.7| ). 



-Cell averages 
- Reconstruction 




-Cell averages 
- Reccnstruction 




-Cell averages 
- Reccnstruction 




Figure 1 : ENO reconstruction of randomly chosen cell averages. 



Remark 1.1. We emphasize that the main Theorem [TT| is valid for any order of ENO reconstruction 
and for any mesh size. It is valid for non-uniform meshes and makes no assumptions on the function 
V, other than that the cell averages v,- must be well-defined, which is guaranteed if e.g. v G L]'p^(]R). 
This is a remarkable rigidity property of the piecewise-polynomial interpolation. 
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Remark 1.2. The stability asserted in Theorem is realized in terms of the reconstructed point- 
values at cell interfaces vf^^^. These are precisely the input for the construction of high-order accurate 
finite volume schemes for nonlinear conservation laws (see Shu ['251), and the relation between these 
values and the cell averages will be the main point of study in this paper. This approach was taken in 
[ITOl . where we use the sign property to construct arbitrarily high-order accurate entropy stable ENO 
schemes for systems of conservation laws. 



Remark 1.3. The proof of both the sign property and the related upper-bound ( |1.8[ ) depends on the 
judicious choice of ENO stencils in Algorithm |1.1[ and it may fail for other choices of ENO-based 
algorithms. In particular, the popular WENO methods, which are based on upwind or central weighted 
ENO stencils, [18, 20, 22J, fail to satisfy the sign property, as can be easily confirmed numerically. 

1.2. ENO interpolation. The ENO algorithm can be formulated as a nonlinear interpolation proce- 
dure. The starting point is a given collection of point- values, {v{xi)}i^z- The purpose of the ENO 
procedure in this context (abbreviated by J^) is to recover a highly accurate approximation of v{x) 
from its point- values vi := v{xi), 

ENO: {v{x,)], ^ Jv{x):=Y,fk{x)\{x), := ^^±^. 

k ^ 

Here, fk{x) are polynomials of degree p—\ which interpolate the given data, 

J^v{Xi) = fi{Xi) = Vi. 

Moreover, the ENO interpolant ^v{x) is essentially non- oscillatory in the sense of recovering v{x) 
to order ff{hP). In particular, since v{x) may experience jump discontinuities, we wish to recover the 
point-values, vT^y^ := //(jc^+i/j) and v+ ■= fi+ii^i+^/i)^ with high-order accuracy, 

v7+U2-v{Xi+l/2-) + vf^,u-v{Xi+i/2 + 



This version of the ENO procedure was used for finite difference approximation of nonlinear conser- 
vation laws in [23|. The ENO interpolant yv{x) is based on the usual divided differences {v[xi, . . . ,;c,+j]},-, 
starting with the grid- values v[x,] = v,- and defined recursively for j > 0. 

Let be the offsets of the ENO stencil associated with grid point x,-. In this case of ENO 

interpolation, these offsets are non-positive integers, corresponding to the integer indices of the pre- 
scribed gridpoints xj. These offsets are selected according to the following ENO selection procedure. 

Algorithm 1.2 (ENO interpolation: selection of ENO stencil). 

Let point values v,_p+i , . . . , v,+p_i be given. 

• Set £i = 0. 

• For each j = 1, — 1, do: 

Uf \v[xi+ej-i,...,Xi+ij+j^i]\<\v[xi+ij,...,Xi+fj+j]\ ^ set £^-+i=£^--l, 
] otherwise i— t- set ij^i=ij. 

• Set fiix) as the interpolant ofv over the stencil {vi+kYkJf^^ ^ ■' 

fi{x) = v[xi+Ej, . . . ,Xi+ej+j] Yl {x-Xi+ij+„) . 

j=0 m=0 

In the following theorem we state the main stability result for this version of the ENO interpolation 
procedure, analogous to the sign property of the ENO reconstruction procedure from cell averages. 
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Theorem 1.2 (The sign property revisited - ENO interpolation). Fix an integer p > \. Given the 
point-values {v,}, let ^v{x) be the p-th order ENO interpolant of these point-values, outlined in 
Algorithm \l.l\ 

= Y,fk{x)\{x), degA(x) <p-l. 

k 

Let v^y2 '^^{^i+^/2~) '^^d =^^('^!+V2+) denote left and right reconstructed point-values 

at the cell interfaces x^j^xj^. Then the following sign property holds at all interfaces: 



(1.9) 




In particular, if the point-values v,+i = v, then the ENO interpolation is continuous across their mid- 
points, v^y2 ~ Moreover, there is a constant Cp, depending only on p and on the mesh-ratio 
—Xj\/\xj —Xj-i \ of neighboring grid cells, such that 

(1.10) 0< ''^+'/^~'''+'/^ <c„. 

V/+1-V/ 

The rest of this paper is devoted to proving the stability of the ENO procedure. We begin with the 
ENO reconstruction procedure in Theorem |1.1[ In Section [2] we prove the sign property ( |1.7[ ) and in 
Section |3] we prove the upper bound ( |1.8[ ). Similar stability results holds for the ENO interpolation 
procedure, dealing with point-values instead of cell averages. In Section|4]we prove the sign property 



for the ENO interpolation stated in the main theorem 1.2 
These results were announced earlier in [|9||. 

2. The sign property for ENO reconstruction 

The aim of this section is to prove the sign property ( |1.7[ ). To this end we derive a novel expression 
of the interface jump, v^y^ ~ ^7+^2' ^ terms which involve {p + l)-th divided differences 

of V, and we show that each summand in this expression has the same sign as v,+i — v,;. 

We recall that at each cell //, the ENO reconstruction is based on a particular stencil of ;? + 1 
consecutive gridpoints, • ■ • ^^i+r+p}, where r = r{i) is the offset of such stencil. 

Notation. We will reserve the indices r and s to denote offsets of ENO reconstruction stencils. We 
recall that these offsets measure the shifts to the left of each stencil, and are indexed at negative 
half-integers, —p + ^/2 < r,s < —1/2, to match the indexing of cell interfaces at half-integers. 

Since this choice of ENO offset depend on the data through the iterative Algorithm |1.1[ we need to 
trace the hierarchy of ENO stencils which ends with the final offset r = r{i). To simplify notations, we 
focus our attention on a typical cell /q, with an initial stencil which consists of the edges at x_i/2 and 
X\i2- The stencil is identified by its leftmost index, ri = —1/2. Next, the stencil is extended, either to 
the left, {jc_3/2, -^-1/25 .^1/2} where r2 = —^/2, or to the right, {x_y2,-<:i/2,.^3/2} with ri = —1/2. In the next 
stage, there are three possible stencils, which are identified by the leftmost offset: rj, = —5/2 corre- 
sponding to {x_5/2, ■ • • ,^1/2}' ^"3 = corresponding to {x_3/2, ■ • ■ ,^^12}^ or ^3 = ~ corresponding 



to {x_i/2, • ■ • ,^5/2}. Stage j of the ENO Algorithm 1.1 involves the stencil of j + 1 consecutive points. 



{xrj,. . . ,Xrj+j}. The series of offsets of this hierarchy of stencils, ri,r2, . . . , r^, forms the signature of 
the ENO algorithm. Note that by our construction, 

r\ = -y2>r2>r2,>--->rp>-p + 1/2, 
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and whenever needed, we set r_i = ro = — 1/2. The stability of ENO will be proved by carefully 
studying such data-dependent signatures. 

The Newton representation of the p-th degree interpolant Fo{x), based on point- values 

V{xrp),V{xrp+i),...,V{xrp+p),is givcu by 

p J-1 

(2.1) Fo{x)=Y,y [^rj , • • • , Xrj+ j] n ~ Xr-j+m) , 

j=0 in=0 

where V[xk,--- ^Xk+j] are the 7-th divided difference of V at the specified gridpoints. Observe that 



in (2.1 ), we took the liberty of summing the contributions of stencils in their "order of appearance' 



rather than the usual sum of stencils from left to right. 

Notation. For notational convenience, we denote the 7-th divided difference of the primitive V as 

D[r,r+j] ■■= V[Xr, • • • ,Xr+j], = ^ _ r = . . . , -3/2, -1/2, 1/2. 

■^r-\- j 



Thus, for example, by ( |1.5[ ) we have £)[-i/2.3/2] = ^ [•'C- 1/2 i-'fi/i 7-^3/2] = (■''^3/2 —^-1/2) (^i ^ ^o)- 



If £>[-i/2,3/2] = 0, or in other words, if vq = vi, then it is easy to see that the ENO procedure will end 
up with identical stencils for Iq and /i, which in turn yields Vy^ = v^,. We may therefore assume that 
-^[-1/2,3/2] 7^ 0, and the sign property will be proved by showing that 



(2.2) Sign property: 



if £>[_i/2,3/2] > then v^^^-Vi^^>0; 
if £>[_i/2,3/2] < then v^2~^i/2-*^- 



To verify (2.2), we examine the ENO reconstruction at cell Iq, given by fo{x) := Fq{x). Differenti- 
ation of (2.1 ) yields 

p j-i ./-I 

/oW = E n i^-^rj+m) ■ 

j=l 1=0 m=0 

The value of /o at the cell interface Xi/2 is then 

p j-i 7-1 p 

(2.3) = /o(x,/2) = £ E n (^V2 - ^0+'") = E ^[rj,rj+j] Yl (■^'/2 " ■^0+'") 

;=1 Z=Om=0 j=l m=0 

The last equality follows from the fact that all but the one term corresponding to / = —rj + 1/2 drop 
out. 

Notation. To simplify notations, we use to denote a product which skips any of its zero factors. 
Thus, for example, a simple shift of indices in (2.3 1 yields Vy, = ^ ^[rj.rj+j] Y\ (^1/2 ^ ^rj+m+i) ■ 

j=l m=— 1 



In an similar fashion, we handle the ENO reconstruction at cell h. Let si, . . . ,Sphe. the signature of 
that cell. Note that rj <Sj + l, since the ENO reconstruction at stage j in cell /i cannot select a stencil 
further to the left than the one used in cell Iq. If rj = sj + I, then the two interpolation stencils are 
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the same, and so Vy^ — Vy^ = 0- Hence, we only need to consider the case Vj < sj. The reconstructed 



value of /i (x) = Fl{x) at x = xi/2 is given by 



j— 1 ra=0 

The jump in the values reconstructed at x = JCi/2 is then given by 



p 



7-1 



J-2 



(2.4) V^2-V,/2= £ ( jT(-^l/2-^.,+m+l) jT. (^1/2 --^O+^+O ) " 



7=1 



m=0 



m=— 1 



The following lemma provides a much needed simplification for the rather intimidating expression 



(2.4), in terms of a key identity, which is interesting in its own right. 



Lemma 2.1. The jump of the reconstructed point-values in (2.4) is given by 



(2.5) 



^1/2-^1/2 



^ ^[r,r+p+l](-^r+p+l " -^r) (-^1/2 " -^r+ra+l ) • 



ra=0 



We postpone the proof of Lemma 2. 1 to the end of this section, and we turn to use it in order to 
conclude the proof of the sign property. To this end, we show that each non-zero summand in ( |2.5| ) 
has the same sign as vi — vq. Since 

(2.6a) sgnj (;C;.+p+i-x,)J^(.JCi/2-x^+,„+i) ) = (-l)''+^"'/^ r = -p + 1/2, . . . , -I/2, 



then in view of (2.2 ), it remains to prove the following. 

Lemma 2.2. Let {ry}^^^ and {sjY^^-^ be the signatures of the ENO stencils associated with cells Iq 
and, respectively, I\. Then the following holds: 

(2.6b) if D[_y2,3/2] > then {-\y+P-'l'D^,,,^p^,^>0, r = rp,...,Sp, 

(2.6c) if Z)[_,/2^3/2] < then {-\Y+P-'l'D^,^,^p^,^<Q, r = rp,...,Sp. 

Since r runs over half-integers, ( 2.6| ) imply that each non-zero term in the sum ( |2.5| ) has the same 
sign as £>[-i/2.3/2]' and Theorem 1.1 follows from the sign property, (2.2). 



Proof. We consider the case (2.6b) where £>[_i/2,3/2] > 0; the case (2.6c) can be argued similarly. The 
result clearly holds for p =1, where = Sp = —^/i. Assuming that it holds for some p>\, namely, 
that {-ly+p-'hD^^^^^p^^^ > for r = rp, . . . , Sp, we will verify that it holds for p-\-\. Indeed, 

(2.7) {-iy+P+y^Dy,^^^p+^^ ^ (_i).+p+i/2 -P[^+L>-+i+P+i]--P[n^-+P+i] 



>0 



for r = Tp, . . . ,5^ — 1, by the induction hypothesis. Thus, it remains to examine £^[^.^+^+2] when 
r = rp+i < rp and r = Sp+i > Sp. 
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(a) The case rp+i = Vp is already included in (2.7 1, so the only other possibility is when r^+i = 
Vp — \. According to the ENO selection principle, this choice of extending the stencil to 
the left occurs when | < |^[rp,rp+p+i] I- Consequently, since (— l)''p+i+/'+V2 — 

(_l)^P+p->/2, and by assumption (-l)'"^+^-V2D[,^ ,^+^+1, = \D 



[rp,rp+p+l] I 



P. 



we have 
>0, 



and ( |2.6b[ ) follows. 

(b) The case Spj\-i — Sp — 1 is already covered in (|2.7|). The only other possibility is there- 
fore Sp^i = Sp. By the ENO selection procedure, this extension to the right occurs when 

l^[.,+i,.,+P+2ll < \Dis,„s,+p+i]\- But (-l)^P+i+P+V2 = _(_i).^.+P-V2 and by assumption, 

|Z)r 



= l^[.>,.>+p+ill'hence 

D[sp + \,Sp+p+2\\ + {-'^y"^''^'^^'D[sp,sp+p+\] 



>0, 



Sp+p- 



-2-Xs, 



and (2.6b) follows. 



□ 



We close this section with the promised proof of Lemma 2.1 



Proof. We proceed in two steps. In the first step, we consider the special case when the two stencils 
that are used by the ENO reconstruction in cells Iq and I\ are only one grid cell apart. Such stencils 
must have the same offset, say rp = Sp = r and in this case. Lemma 2.1 claims that Vy^ — equals 



P~i 



(2.8) 



m=0 



Indeed, the interpolant of V{xr+i),. ■ ■ ,V{xr+p),V{xr), assembled in the specified order from left to 
right and then adding V{xr) at the end, is given by 

p-i j-i p-i 

Similarly, the interpolant of V{xr^i),. . . ,V{xr^p^\), assembled in the specified order from left to 
right, is given by 

p-\ j-l p-i 

-^1 W = ^ ^[r+l,r+;+I] ~ ^r+l+m) + ^[r+I,r+p+I] ~ ^r+l+m) 

j—0 m=0 m=0 



(cf. (2.1 )). Thus, their difference amounts to 



p-i 



Fi{x) -Fo{x) = (£>[r+Lr+p+l] -D[r^r+p]) Y[ i^^^r+l+m] 

m=0 

p-l 

= ^[r,r+p+l] {Xr+p+l — ^r) Y[ — Xr+l+m) ■ 



m=0 
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which reflects the fact that Fq and Fi coincide at the p points Jc^+i, . . . ,Xr-\-p. Differentiation yields 

p-ip-i 

/l W -/oW = D[r,r+p+l] {xr+p+1 ' ^r) Yl i^-^r+l+m) ■ 

1=0 m=0 

At X = xi/2, all product terms on the right vanish except for I = —r— 1/2, since Zi/2 belongs to both 
stencils. We end up with 

p-l 

m=0 



This shows that (2.8 1 holds, verifying Lemma 2.1 in the case that the stencils associated with Iq and 
/i are separated by just one cell. 

In step two, we extend this result for arbitrary stencils, where Iq and /i are associated with arbitrary 
offsets, rp < Sp. Denote by F^j^ the interpolant at points xj, . . .,Xj+p, so that Fq = F^''''^ and Fi = 
p{sp+i} _ Using the representation from the first step for the difference of one-cell shifted stencils, 
(^^■{'■+1} {xi/2), we can write the jump at the cell interface as a telescoping sum. 



(/i -/o)M = (/^^''+^^ -/^^'^^) M 



p-i 



^ ^[r,r+p+l]i^r+p+l — Xr) ^ {xi/2 — Xr+m+\) 



m=0 



and ([23]) follows. □ 
3. The relative jumps in ENO reconstruction are bounded 



In this section, we will prove ( |1.8[ ), which establishes an upper bound on the size of the jump in 
reconstructed values in terms of the jump in the underlying cell averages. We need the following 
lemma. 

Lemma 3.1. Let rp,Sp be the (half-integer) offsets of the ENO stencils associated with cell Iq and, 
respectively, Ii. Then 

(3.1a) ^^vlZ!±ll(_i)^+.->/2 < . = rp^ ...^Sp, 

^[-1/2,3/2] 

where the constants Cr,p are defined recursively, starting with C^.i = 1, and 

2 , . 

(3.1b) = — max(Qp,Q+i,p) V r. 

Xr+p-\-2 Xy 



The quantity on the left in p.la[ ) was shown to be bounded from below by zero in Lemma 2.2 , here 
we prove an upper bound. The constants Cr,p only depend on the grid sizes \Ij\. 

Proof. The result clearly holds for ;? = 1 . We prove the general induction step passing from p\—^p+\. 
Using the recursion relation 

_ D[r+l,r+p+2\ -D[r,r+p+l] 
^[r,r+p+2\ — ; ^ •— Xr+p+2 —Xr, 
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we have 



(3.2) 



~ ^[-1/2,3/2] ^ \ ^[-1/2,3/2] ^[-1/2,3/2] / 



< 



Cr+l,p~\'Cr,p 

Ax 



r — rp, . . . ,Sp 1. 



We turn to the remaining cases. 

(a) As in Lemma 2.2 if r^+i = then the induction step is already covered in (3.2), so the 
only remaining case is r = r^+i = — 1, corresponding to Lemma 2.2 a). In this case, 

\D[r,r+p+i] \ < l^[r+i,r+p+2]l'hence 

D[r,r+p+2] ^ ^y+p+l/2 ^ 1 D[r+l,r+p+2]-D[r,r+p+l] ^ ^y+P+l/l 



D 



-I Ml] 



Ax 



/2,y2] 



< 



2 D[r+i,r+p+2] I + +1/2 2Ch-1j. 

Ax £>[_i/2 3/2] Ax: 



(b) If =Sp — \ then the induction step is already covered in (3.2 1, so the only remaining case 
is r = Spj^x = Sp, corresponding to Lemma 2.2 ;b). In this case, \J^\r^\^r^p^2"\ \ ^ l^[r,r+p+i] I' 
hence 



-1/2,3/2] 



< 



1 D[ 


r+l,r+p- 


Ax 


D 


2Z)[ 




Ac D 


-'/2,3/2] 



1/2,3/2] 



-(^_1)'"+P+V2 



(-1) 



2C, 



□ 



Using the explicit form ( 2.5) of the jump v^^ — Vj^^, we get the following explicit expression of the 
upper-bound asserted in ( |1.8| ). 

Theorem 3.2. Lef v^, anJ Vj^^ point-values reconstructed by the ENO algorithm ^. \^ at the cell 
interface x = X1/2+ and, respectively, x = Xxj^—- Then 



^^/22V 

Vl -vo 



1 



-1/2 



^3/2-^-1/2 ,= _p+I/: 



P-1 



(-^r+p+l ~-''^f)^p^ (.''^1/2 — 



ra=0 



Proof. Let rp,5p be the offsets of the ENO stencils associated with cell /q and, respectively, l\. By 



Lemmas 2.1 and 3.1 we have 



^1/2-^1/2 



1 ^ T> 



Vl-Vo -^3/2-^-1/2 r= 



[r/+p+l] ^ yy^p-Xjl 



MA 



p-1 



m=0 



< 



1 



p-1 



< 



X3/2-JC_l/2 r=,p 
1 

X3/2 X_i/2 r=-p+l/2 



(-"^r+p+l ~-''^) ('"'^1/2 ^■''^r+m+l) 



ra=0 



p-1 



(■"^r+p+l —Xr) (xi/2 — .'C^+m+l) 



m=0 



□ 
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When the mesh is uniform, |/, | = h, the expression for the upper bound C can be calculated explic- 
itly. The recursion relation p.lb[ ) yields Cr^p = ^ ^ ^ ^ — TTVT' ^^'^ coefficient of the {p+ l)-th 



order divided differences in (2.5 ) is 



p-i 



?n—0 



hP{p+l){-r-l/2)\{p + r-y2)\. 



Thus, we arrive at the following bound on the jump in reconstructed values. 

Corollary 3.3. Let v^, and v^^ be the pointvalues reconstructed at the cell interface x = X\j2+ and, 
respectively, x = Xxj^—, by the p-th order ENO Algorithm |1.1| based on equi-spaced cell averages 
{vk}h Then 



(3.3) 



^V2-^l/2 

Vi -vo 



<Cp = 2P-'-'Y,k\{p-k-l)\. 



k=0 



Table 1 shows the upper bound p.3[ ) on (vy^ — Vy^) / (vi — vq) for some values of p. Returning to 
Figure T we see that although the jumps at the cell interfaces can get large, they cannot exceed Cp 
times the size of the jump in cell averages, regardless of the values in neighboring cell averages. 



p 


Upper bound Cp 


1 


1 


2 


2 


3 


10/3 = 3.333... 


4 


16/3 = 5.333... 


5 


128/15 = 8.533... 


6 


208/15 = 13.866... 



Table 1. 



It can be shown that the bound Cp given in Theorem 3.2 is sharp. Indeed, the worst-case scenarios 
for orders of accuracy p = 2,3,4,5 are shown in Figure [2| The mesh in this figure is x,_|_i/2 = i, and 
the cell averages are chosen as 



V,: 




if i is odd 

if / is even and i < 4 
if i is even and i > 4. 



The number 10^^*^ is chosen at random; any small perturbation will give the same effect. This pertur- 
bation ensures that cells /, for i < 4 interpolate over a stencil to the left of the cell interface x = X/^j^xj^, 
and cells with / > 4 to the right of it. We see that for each p, the jump at x = 4 is precisely the bound 
Cp as given in Table [T| 



4. ENO INTERPOLATION 



The stability proof of ENO interpolation stated in Theorem 1.2 can be argued along the lines of 
those argued in Sections [2] and [3| We therefore only sketch the arguments without details. 
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p = 2 



p = 3 




« / 


/ / 


Cell averages 

Reconstruction 

} o 









9 1 



Figure 2: Worst case cell interface jumps for p = 2,3,4,5. 



4.1. The sign property for ENO interpolant. We focus on the jump across the interface at .x; = jci/j. 
Let and {njl^^j be the ENO stencils associated with gridpoints xq and, respectively, xi. 

Recall that in this case of ENO interpolation, these offsets are non-positive integers, —p + 1 < ij.nj < 
0. Our first key step is to compute the jump at the interface point Xi/2 in the case where the ENO 
procedure are separated by just one point, namely, £p = Up. We let abbreviate the divided 

differences v[x/, . . . As before, we denote the reconstructed values at the interface x = X\j2 by 



1/2 



Jv{x\j2-) and v^i^ = yv{xi/2+)- 



Lemma 4.1. !/£„ = np = I for some i G — Nq, then 



m=\ 



"1/2 "1/2 

By assembling a telescoping sum of several such stencils we obtain 
Corollary 4.2. For general £p < np, we have 



(4.1) 



12 d[e/+p+i]ixe+p+i 



p-i 

m=l 



Since vi — vq = {xi —xo)d[Q ^, we wish to show that the jump in reconstructed values at the cell 
interface has the same sign as ^i?[o,i]- To this end we show that each summand in ( |4.1[ ) has the same 
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sign as d^Q iy Indeed, since 



p-i 



sgn (xe+p-Xi) Y[(^y2-X£+m) = (-1) 



-p + l<£<0, 



it suffices to prove the following: 

Lemma 4.3. If I p, rip are selected according to the ENO stencil selection procedure, then 

if J[o,i]>0 then (-l)^+^+^rf[,,,+p+i]>0; 
if J[o,i]<0 then (-l/+^+^,,,+p+l]<0, 

4.2. Upper bounds on the relative jumps for ENO interpolant. Next, we show the corresponding 
upper bound on v^^ — Vy, for ENO reconstruction with point values. 

Lemma 4.4. Iflp, np are selected according to the ENO stencil selection procedure, then 



0< 



[0,1] 



where c^n are defined recursively, starting with Cy.i = 1, and 



<^e,P+i 



p , . . . , np , 



■max(c^,p,c^+i,p). 



p 


Upper bound Cp 


1 


1 


2 


2 


3 


3.5 


4 


6 


5 


10.375 


6 


18.25 



Table 2. 



For simplicity we assume that the mesh is uniform with mesh width x^+i —xj = h. It is straight- 
forward to show that ce„ = (2//i)^^^ ^/p- Moreover, the coefficient of the {p + l)-th order divided 



differences in (4.1 ) is 



p-i 

(Xf+p -Xl) W (Xi/2 



hPp 



p-l 
m=l 



m 



Thus, we arrive at the following bound on the jump in reconstructed values. 

Theorem 4.5. Let ip,np be selected according to the ENO stencil selection procedure, and assume 
that the mesh is uniform. Then 

,+ 



^1/2-^1/2 

Vi -vo 



< c 



p 



2P' 



-1 



1 



(P 



-IV ^ 



p-l 

m=l 



2 — i 



-m] 



Tablej2jshows the upper bound on (vy^ — Vy^) / (vi — vq) for p < 6. As for the ENO reconstruction 
procedure in Section [3| it may be shown that these bounds are sharp. 
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5. Conclusions 

We show that the ENO reconstruction procedure (from cell averages or point values) is stable via 
the sign property, namely the jump in reconstructed values at each cell interface have the same sign 
as the jump in the underlying cell averages (point values). Furthermore, we obtain an upper bound 
on the size of the jump of reconstructed values in terms of the underlying cell averages (point values) 
at each interface. Both results hold for any mesh {x/^i/j}/- In particular, the results hold for non- 
uniform meshes. In addition, both results hold for any order of the reconstruction, i.e, any degree for 
the polynomial interpolation. No extra regularity assumptions on the underlying Lj^^ function v are 
needed. 

The proof of both the sign property and the upper jump bound depended heavily on the formula 



( |2.5[ ), which gives the cell interface jump in terms of r^, Sp and the {p + l)-th divided differences 



of V. This formula is completely independent of the ENO stencil selection procedure, and hence 



holds for all interpolation stencils. On the other hand. Lemma 2.2 (and Lemma 3.1 for the upper 
bound) is a direct consequence of the ENO stencil procedure. Therefore, we cannot expect that other 
reconstruction methods satisfies a similar sign property. In particular, the WENO method, using the 
stencil weights proposed in [|T8l l22| . will in general not satisfy such a property, a fact that is easily 
confirmed numerically. This leaves open the question of the existence of stencil weights that make 
the method satisfy the sign property. Of the second-order TVD reconstruction methods (see [|26l ). 
only the minmod limiter satisfies the sign property. 

The stability estimates presented in this paper do not suffice to conclude that the ENO reconstruc- 
tion procedure is total variation bounded (TVB). In particular, the jump in the interior of a cell can 
be large. However, the sign property enables us to construct arbitrarily high-order entropy stable 
schemes for any system of conservation laws. Furthermore, the sign property together with the upper 
bound allow us to prove that these entropy stable scheme converge for hnear equations. Both results 
are announced in ^ and presented in a forthcoming paper [[TOll . 
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